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. The rate constants for recombination and exchange reactions are calculated using the flux corre- 

' lation approach with a general form of the Boltzmannized flux operator, which can simultaneously 

' describe the Kubo and traditional half-split iorms. First, we consider an exactly solvable model, i.e., 

^Nj , the free particle case, in terms of a new scaling function. Next, as a non-trivial case, we study the 

recombination and exchange reactions at very high pressure. Since the rate constant is calculated 
^ l' by Laplace transform of the flux correlation function, the result depends on how the Boltzmannized 

, flux operator is chosen. We flnd that a choice of the flux operator affects the rate constant consid- 

■ erably. For the recombination reaction, the ratio of the rate constant in the half-split form to that 

CO ' in the Kubo form approaches zero in the high pressure limit. 
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I. INTRODUCTION 



To study a chemical reaction at the most detailed level, it is necessary to treat the action of the time evolution 
operator, exp{—iHt) {H is the Hamiltonian of the system), ^ onto a (given) initial wave function or density matrix, 
^ which is a function of total energy E and total angular momentum J. Such quantum scattering calculations have 
actually been studied for simple chemical reactions.^ A time-dependent scattering formalism based on the 5- matrix 
^ Kohn variational approach^ or a coupled channel method in hyperspherical coordinatesi^ has usually been used to 
Q , evaluate the quantum reactive scattering cross sections. The number of open vibration-rotation channels, however, 
■ increases very rapidly as thermally accessible collision energy becomes high. Therefore, at high energy the exact 
quantum state-to-state calculations would not be feasible even for a simple reaction. Furthermore, if it is only the 
rate constant that is desired, such a complete calculation for all statc-to-state information is not economical. 
CIh a totally different approach for treating chemical reactions is that based on a correlation between quantum flux 
operators in the statistical thermodynamics. In the early 60's, Yamamotoi first formulated an exact expression 
for the rate constant as an application of the general statistical mechanical theory of irreversible process, which 
J> was established by Kubo et ali^ and Moriii (See also RefH.) Independently, Miller et ali2*i2iii have developed a 
00 method for the rate constant using a time integral of the flux-flux autocorrelation function, which is also exact in 
OO \ the limit that the dynamics is extended to < cxd. The feasibility of this approach depends on how to evaluate the 
Op • time evolution operator for the system, and the correlation is usually calculated using the Feynman path integral 
\ techniqueiiSii^ Although this method is powerful and convenient to obtain the rate constant, the direct application to 
, large, complicated systems is still difficult because of the notorious sign problem>i^ii^ One approach to improve this 
( <~>) ' situation is a filtering or smoothing technique^^ such as the stationary-phase Monte Carloiii Recently, semiclassical 
(SC) approaches, implemented via the initial- value representation (SC-IVR), have received a rebirth of interest. (Van 
Vleck^'^ first discussed a drastic approximation based on the semiclassical picture.) A number of studiesi^iSSiSi have 
^ , been carried out along SC-IVR and have demonstrated the capability of these approaches to describe various quantum 
effects in chemical reactions. 

' It has also been shown how a quantum mechanical version of the Lindemann mechanism for coUisional rccombina- 
Qh' tion, where the process is affected by the bath gas M, 



A + B ^ AB*, (1) 

AB* + M AB + M, (2) 

can be handled by the flux-flux autocorrelation function for the A ~ B collisioniSSiSil'* Some applications of this 
theory are listed in Refs l25ll2^ It is furthermore possible to generalize the formalism to include chemical reactions as 
well as recombination: 

A + BC ^ ABC* AB + C, (3) 

ABC*+M ABC + M. (4) 



t We use the natural unit, i.e., ti = c = 1. 
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Equations Q and Q simultaneously describe the recombination process {A+BC — > ABC) and the exchange reaction 
{A + BC — + AB + C).^^ This method has been applied to the interesting (combustion) reactions (O + OH ^ H + O2) 
and the recombination reactions (O + OH + M HO2 + ^/ ^ H + O2 + M) , which are very important in atmospheric 
chemistryi^iS^ 

As mentioned above, there are two different ways to calculate the rate constant for chemical reactions, i.e., the 
Yamamoto approach, in which the Kubo form of the flux operator is used, and the Miller approach, in which the 
traditional half-split form is chosen. The two approaches can provide the same result to the rate constant because 
it is given in terms of the integral of the flux-flux correlation function with respect to timeAii However, the shapes 
of the correlation functions calculated by the two approaches are quite different from each other. Therefore, for 
the recombination and exchange reactions (like Eqs.© and the rate constants in the two approaches could be 
different, because they are given by Laplace transforms of the flux-flux correlation functionsi^ The purpose of this 
paper is to demonstrate the effect of different choices of the Boltzmannized flux operator on the rate constants for 
recombination and exchange reactions. 

First, we summarize the correlation function method briefly in Sec|n] The expressions of the correlation functions 
in the Yamamoto and Miller approaches can be unifled using a general form of the Boltzmannized flux operator. The 
difference between the two approaches is discussed explicitly. The rate constants for recombination and exchange 
processes are also studied. In Sec lIIII we consider an exactly solvable model, i.e., the free particle case, in terms of 
a new scaling function. In Scc lIVI as a non-trivial case, we study the recombination and exchange reactions at very 
high pressure. Finally, the summary and conclusion are given in SeclVl 



II. FLUX-FLUX CORRELATION APPROACH TO RATE CONSTANTS 

The quantum mechanically exact expression for a thermal rate constant k{T) can be written in terms of the flux 
correlation functional-'- 

fc(T) = go(T)-Mim a(r;t), (5) 
where Qo{T) is the reactant partition function per unit volume and Cs(r; t) is the flux-side correlation function defined 

by 

Csir;t)^tr[e~^^F{r)'P]. (6) 

Here, F{r) is the bare flux operator given by 

F{r)^i[H,h{s)l^r, (7) 

with h the Heaviside step function and s the reaction coordinate - see Fig^ h{s) takes the value of 0(1) in the 
reactant (product) side of the dividing point r. In Eq.®, V is the projection operator deflned by 

r = e"'h{p)e-'"\ (8) 

with p the momentum operator. 

For the rate constant, two different approaches have been proposed so far: one is the Yamamoto approach)4 which 
is based on the linear response theory (or the so-called Kubo formula))^ and the other is the flux-flux autocorrelation 
function method developed by Miller et alAiSiii In the Miller approach, Eq.® are modifled by following two useful 
facts: the first one is made by noting that V and H commute each other (i.e., [P, H] = 0) in the limit t 00, and the 
second is to replace h{p) by h{s) in Eq.®, which is also correct in the limit t — > 00. By virtue of these modifications, 
the flux-side correlation function in the Miller approach can be rewritten as 

Cf{r-t) = iY[F^\(3,r)e'"'h{f)e-'"% (9) 

where the superscript M stands for "Miller". Here, F*^(/3,r) is the half-split Boltzmannized flux operator 

F^i (/3, r) = e-P"/^F{r)e-^"/\ (10) 



t For simplicity, we consider a reaction in one-dimension. 
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Using Eq.lO, the rate constant can also be expressed in terms of the flux-flux correlation function C{r, r; t) as 

/•oo 

Qo{T)k{T)^ / dtC^'ir,r;t), 
Jo 

where 

C^(r, r; t) = tr[F^^{p, r)e'"' F{r)e-'"% 
On the other hand, the linear response theory gives the rate constant 

/•OC 

Qo{T)k{T)= / dtC''{r,r;t), 

where 



(11) 



(12) 



(13) 



(r, r; t) = tr[f^>^ (/3, r)e*"'F(r)e 



(14) 



with 



P Jo 

Here, the superscript Y stands for " Yamamoto" . 

If a general form of the Boltzmannizcd flux operator is introduced bjjSl 

1 /•(1 + k)/3/2 

F^iP, r) = 4 / dX e-(^-^)^/'(r)e-^^, 

J(1-k)/3/2 



(15) 



(16) 



where k is a parameter (0 < n < 1), the correlation functions in the Miller and Yamamoto approaches can be unified 



as 



(17) 



It is easy to check that in the limit k 0(1) Ea. (|17|l reproduces the correlation function in the Miller (Yamamoto) 
approach. 

Using Eq. jT)) , the integral with respect to A in the general form of the Boltzmannized flux operator can be performed: 



F-iP,r) 



Kf3 



(18) 



Combining the partial Boltzmann operator, e (i±«)/5^''/2^ with the real-time evolution operator e^'^*, we obtain two 
complex-time evolution operators: e~*^*+ and e~*^*" with t± — t — i{l zL k)(3/2. Then, the correlation function is 
rewritten as 



C^(r,r;t) = —3 tT\F(r)e'"''-^h(f)e-'"'- 



(19) 



where 9 stands for taking the imaginary part. Performing the trace operation in Ea. (|19|l . the correlation function 
reads 



C''(r,r';t) = tr\F(r)e'^*+h(?)e-'"'-] 



1 



t(s\e~'"'+\s')*^(s\e-'"'-\s') 



ds' 



^{s\e-'"'~W}^{s\e-^"'^\sr 



(20) 



where (s|e *^*|s') is the propagator in the coordinate representation. The correlation function C'^{r^r]t) is given in 
the limit r' — > r. 
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It is possible to generalize the flux correlation approach to treat recombination and exchange reactions. ^iS422iS4i2Si2& 
Let us consider the reaction of ^ + BC AB + C and ABC in one-dimension (see Eqs.Q and 10} and Fig|2Il. The 
rate constants for the exchange {A + BC AB + C) and recombination {A + BC ABC) reactions are again given 
by the average of the flux F{r) and the projection operator V over the Boltzmann distribution, where F{r) describes 
the flux at the reactant dividing point r (see Fig|2l. Similarly, we can define the flux operator F{p) at the product 
dividing point p. 

Because the probability of the system experiencing a deactivating {ABC* + M ^ ABC + M) collision with the bath 
gas M can be evaluated by 1 — e"''* at time t {rj describes the frequency of deactivating collisions and it depends on 
pressure P and T of the bath gas), the recombination probability is proportional to 1 — e^''"^, where r is the time the 
trajectory (it is on r at i = 0) is in the compound region (see Figl^J- On the other hand, for the exchange reaction, 
the probability is given by e~'^'^^ , where is the time the trajectory exists through the point p. 

Using these probabilities and the correlation function, the rate constants for the recombination and exchange 
reactions are finally given byS^ 

Qo{T)k^,,{T,P) = / dte-^'[C''{r,r;t)-C''{r,p;t)], (21) 
Jo 

/>oo 

Qo{T)K.c{T.P) ^ dte-^'C'^{r,p;t), (22) 

where the relation 

roo poo 

dtCir^r-t)^ dtC''{r,p;t) (23) 

holds. This ensures that the recombination rate vanishes in the limit 77 — )■ 0. 

Because the rate constants for recombination and exchange processes are calculated by the Laplace transforms 
of the correlation functions, it is clear that the evaluated results depend on k, that is, the shape of the correlation 
function affects the rate constantsi^S Of course, the integral of the correlation function with respect to time itself does 
not depend on k. 



III. SOLVABLE MODEL - FREE PARTICLE CASE 



In this section and the next section, we calculate the rate constants for recombination and exchange reactions using 
the Feynman path integral technique]— ii^ A huge calculation is usually required to obtain the exact matrix elements of 
propagators for a realistic system. Furthermore, it is necessary to take some approximations and numerical techniques 
to perform iti^7ii?i29i2ii?9i?ii?2 Because the aim of this paper is to show how the rate constant for recombination or 
exchange reaction depends on k (i.e., a choice of the Boltzmannized flux operator), it would be more intuitive and 
useful to flrst consider a simple system rather than a complicated case. We here study an analytically solvable model, 
i.e., the free particle case in one-dimension, and discuss a non-trivial case in the next section. 

If no potential acts on the system (see FiglSj), the propagator in the coordinate representation is easily calculated 
by the path integral«i^ The flux-flux correlation function Ea. H2(J|l then gives 



where D±=t^ + {l± nyP^/A, X 



A = 



B = 



1 



exp 



X [{2tA- K[3B)iiinX 

Kmtp'^d'^ I2D+D^ and 
1/2 



m/3d2(t2 + (i„K2)/32/4)' 



{2tB 



k(3A) cos X] 




(24) 



(25) 
(26) 



with d the distance between r and p and m the reduced mass. Note that, as it should be, the correlation depends on 
only the distance d and is independent of the positions r and p. 
If K = 1 (the Yamamoto approach) , Eq. (|24|) gives 

m(3d^ 



Q{r,p;t) 



1 



271(3^ y/2t(W+W) 



exp 



2(i2 + /32) 



(Vt^T^ + i)^/' sinX' + (Vt' + /32 - t)3/2 cosX' 



(27) 
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where X' = m[3^(P/2t{e + (3^), and 



C(r,r;i) 



27r/32V2t(t2+/32)i/2- 

On the other hand, taking the hmit k — > in Ea. H24|) . we obtain the correlation functions in the Miher approach 



(28) 



and 



47r(t2 + /32/4)3/2 



13 2rat^(P 

2 ^ l^TWJ^. 



exp 



2(i2 + /32/4) 



B7r(t2+/32/4)3/2- 



(29) 



(30) 



Now it is very convenient to introduce a new scahng function 

1 



: exp 



a;2 + (l-K2)/4) 



D+D^ 



2Ky/D+D 

X [(2xA- kB) sin X + {2xB + kA) COS X] , 



(31) 



where x{— tj (3) and c{— m<P 12(3) are dimensionless variables. Since the scaling function itself is a function of x and 
c, it is also dimensionless. Here, D± = a;2 + (1 ± k)^/4, X = kcx/D+D-, and A and B are given by replacing D± 
with Z)± and t ^ x in Eas. lf^ and respectively. If we set r = p (or d = 0), we obtain 



S^^{x) = 27rP^CS{r,r;t) ^ 



1 



2k^D+D^ 



-.{2xB + kA), 



(32) 



which is a function of only x because c = 0. 

We show these scaling functions in Figs 01 and [S] Note that the integral of the scaling function with respect to 
X does not, of course, depend on k. In Fig^jthe scaling function for k — 1 (Yamamoto approach) is divergent like 
~ I/a/S at a: = 0. Taking the limit a; — > implies that t approaches but f3 is finite or T is very low (/3 — > oo) but t 
is finite. In the Miller approach (k — 0), the scaling function is finite at x = 0: S'^^°(0) — 2. For arbitrary k, we find 
S'^^(O) — 2/\/l — k2. By contrast, in the limit x — > oo (which corresponds to the case where T is very high {(3 0) 
but t is finite or t — s- oo but /3 is finite), the scaling function does not depend on k and behaves like l/4x^. 

In Fig[Slwe take c — 1.0 to illustrate the scaling function S^p clearly. It is remarkable that the interference effect in 
the correlation is fully taken into account in the Yamamoto approach (k — 1). Hence the scaling function oscillates 
very rapidly near the origin. Contrastingly, in the case of k ^ 1 the interference is averaged and the shape is very 
smooth. The value of the scaling function at x = is given by (2/Vl - K2)exp(-4c/(l - k^)) for k ^ 1. The 
partial-split form of the Boltzmannized flux operator (i.e., < k < 1) is thus very convenient for actual computations. 
For large x the scaling function is again independent of k and behaves like '--^ (1 + 8c) /4a;'^. 

Next, we study the rate constants for recombination and exchange reactions, which are given by Eas. (|21|l and H22|l . 
respectively. Using the scaling functions, we first calculate the (rr)- and (rp)-rate constants, which are, respectively, 
defined by 



Qo{T)K 



1 



2^/? Jo 

QoiT)Kp = dxe-'^^S^pix.c), 



(33) 
(34) 



with a{= (3ri) a dimensionless variable. Then, the rate constants are given as k^^^ — ^rp ^exc — ^rp' 

In order to convert the collision frequency rj to more familiar variables, we approximate the coUisional deactivation 
rate constant by an expression given by the hard sphere collision theory. Furthermore, if one uses the ideal gas 
expansion, the frequency can be expressed by2^ 



V = kdeact[M] ^ P 



2000 



T 



X 10" 



(35) 
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with ?7 in fs ^, P in Pa and T in K. Then, we find 

~ /19 ^/ 

T3/2 

The factor c is also converted as 



a = prj^ 3.42 x x IQ-^. (36) 



c = — = 0.0103 X A^Td\ (37) 
2p 

with Am the reduced mass of the system in atomic mass units and d in A. 

Now we are in a position to show our results for the free particle case. We define the ratio of the rate constant for 
arbitrary k to that for k = 1 (Yamamoto approach): 



R'^^ia) = KJK=\ (38) 



R^,Ja,c)^KJK-\ (39) 



The ratio is a function of only a, through which it depends on T and P, while the ratio i?^^ depends on a and c. 
In FigsEland[3 we illustrate the two ratios for the free particle case. Here, we choose < a < 1, which, for example, 
covers the region of P < 1.5 GPa at T ~ 300 K. In FigO c is fixed to be 20 because the ratio is not sensitive to it in 
the region of < a < 1. In i?^^, the ratio decreases linearly and it is about 0.9 (10% reduction) at a = 1 in the case 
of K = 0. Note that a large value of a corresponds to the case of high P and low T, and that the ratio at large a is 
mainly determined by the correlation function (or the scaling function) at a; ~ (see Fig^J. By contrast, in i?"^ the 
deviation of the ratio from unity is not large, and the ratio for k = is enhanced by about 4 — 5% at a = 1. 

Combining the (rr)- and (rp)-rate constants, we can calculate the ratio of the recombination rate constants, -RJ^e^ = 
krecl^rTc ■ Bccausc the (rp)-rate constant itself is, however, much smaller than the (rr)-rate constant, the ratio R^ec 
is very close to i?^^. We find that is again about 0.9 at a = 1. On the other hand, the ratio for the exchange 

process is given by i?^^^ = K^c/K7c = Kp- 

IV. RATE CONSTANTS AT HIGH PRESSURE LIMIT 

As a non-trivial example, we consider the rate constants for recombination and exchange reactions at very high 
pressure. As we have seen in the previous section, it is very notable that at short time (or small x) the correlation 
function (or the scaling function) in the Yamamoto approach is quite different from the other case (k ^ 1). The 
rate constant is calculated by the Laplace transform of the correlation function (sec Eqs.|(2U ^"^^ 122|l ). in which the 
integrand involves the factor exp(— yyi) and 7] is proportional to pressure. Thus, at sufficiently high pressure the rate 
constant may be determined by the correlation function at short time (t < 77^^). 

For a short time, the paths which must be considered in evaluating a potential V never move very far from the 
initial position. Thus, to a first approximation, we can expand the potential V around the average of the initial (s) 
and final (s') positions: 

V{x{t)) ~ V{x) + [x{t) - x]V'{x) + - xYV"{x), (40) 

where x = {s + s')/2. Here we ignore terms of higher order than the second derivative of V . By virtue of the quadratic 
form of x{t), we can easily find the propagato r ^^i^^ 



2-1X1 sin T 



1/2 



itV{x) iS 



e 'e 



(41) 



where 



(l+eosr)(.-.r- ZS_ ^ (42) 



4 sin T 2mu^ sin r ' 

with r = utt, f —V'{x), uP' — V"{x)/m and F{t) — 2(1 — cost) — rsinr. Because the potential is expanded up to 
0{V"), Eq.jni) is valid up to 0{p) or 0{uj'^). 
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To take the temperature effect into account, it is again necessary to introduce complex-times. Because the correlation 
function is given by (see Eas. ((TC|l and ((T7|l ) 

1 /•(1+k)/3/2 . ^ - - „ 

C'^[r,p-t) ^ — dXtT[e-^^-^^"F{r)e-^"e'"'F{p)e-'"% (43) 

'^P J(l-K.)B/2 



1 
1 



(1-k)/3/2 
K/3/2 

dA tr[i^(r)e^^(*+^('3/2+A))^(p)g-zff(t-^(/3/2-A))]^ ^44) 

-K/3/2 
K/3/2 

dA C"'(r,p,A;t), (45) 

-K/3/2 



we define new complex-times as t±(A) — t zt i{(3/2 zL X). Because of very high pressure (or short complex-time), it is 
clear that the (rp)-rate constant becomes very small compared with the (rr)-rate constant. Therefore, we here focus 
on the (rr)-rate constant and do not consider the (rp)-rate constant. Then, using the propagator Ea. (|41|l we obtain 



C^ir,r,X;t) 



1 



4^[i+(A)<_(A)]3/2 



P 0' 



•3 



2 24to 



V"{r) 



e-'^^W, (46) 



up to 0{V"). The correlation C^^r^r^X^t) is completely determined by V and V" , while the first derivative V gives 
a correction to the correlation for the (rp)-rate constant. Because the quantum correction due to the potential does 
not depend on A explicitly (sec Ea. H46l) '). the correlation function C"'(r, r;t) can be described in terms of that for the 
free particle case. It is then written best in the form 



C''{r,r;t) = exp 



•2 



xCS{r,r-t). (47) 



In a typical chemical reaction, the interatomic distance and forces usually range over one or two angstroms, and 
hence one can expect that the change of the potential is small while the motion of the order of lA has been achieved. 
Because, in the ideal gas model, the root-mean-square speed of a molecule is usually estimated as v ^ 1.6-\/r/A,„ x 10^ 
(m/s), the particle moves by as much as Ax ~ 3.5{T/Py/Am) x 10"'' (m) for the short period of rj^^. Thus, if Ax < 1 
A, the expansion of the potential Ea. H4U|) is justified and we can use the correlation function Ea. H47|) . If we assume a 
typical mass of 20 and T - 300 K, we find the condition of P > 0.24 GPa. 

Since the correction due to the potential in Ea. H47() is common to the correlation function for arbitrary k, the ratio 
of the (r'r)-rate constant for arbitrary k to that in the Yamamoto approach (k = 1) is again given by the ratio for 
the free particle case: i?^^(a) = k'^^./k'^^^. (Note that R'^^^ ~ iJJ!^ because at high pressure the (rp)-rate constant is 
negligibly small.) In Fig|Sl we show the ratio in the region of 1 < a < 10 (cf. Fig|SJ). At large a the ratio is 
reduced considerably, which implies that at very high pressure the (rr)-rate constant for k ^ 1 is much smaller than 
that in the Yamamoto approach. 

At high P, the rate constant for k 7^ 1 may be estimated by 

Q^{T)Kt ^ - X ^;/i(0) / dx e-"- = (48) 

where Z is a (dimensionless) constant which involves the correction due to the potential. By contrast, in the case of 
K = 1 the rate constant may be given by 

, _ Z . ^„._, f . e 



{T)k^=' ~ — X [V^S:=\x)],^o I dx^^ (49) 







The ratio at high pressure is thus given as 



Vl - k2 V 7ra 



which approaches zero in the limit P — > 00. At a = 10 (which, for example, corresponds to the case of T ~ 300 K and 
P - 15 GPa), Eq.(|Sni gives = 0.5(0. 58)[0. 76] for k = 0(0.5)[0.75]. The value for k = very well agrees with the 
numerical result (see Fig|Sl), while the value for k ~ 0.75 is a little larger than the numerical one. The approximation 
Ea. (|48|l may not work well for large k because the peak of the scaling function at a; = becomes sharp near k = \ 
(see Fig01). In the present calculation, we find that the (rr)-rate constant in the Yamamoto approach is considerably 
different from the result calculated by the partial-split form. At a ^ 10 the rate constant for k = 1 is about twice as 
large as that for k = 0. 
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V. SUMMARY AND CONCLUSION 



The exact quantum mechanical expression for thermal reaction rates can be formulated by the linear response 
theoryj^iSii which Yamamoto first discussed in the early 60's.'*> Later, Miller et al&iSiii have independently proposed 
a more convenient way, i.e., the flux-flux autocorrelation function method, to perform numerical computations. The 
Miller approach can provide the exact rate constant in the limit that the dynamics of the system is extended to 
t —^ oo. Using a general form of the Boltzmannized flux operator, we have unified the two approaches and studied the 
rate constants for thermal exchange and recombination reactions. Because they are calculated by Laplace transforms 
of the flux-flux correlation functions, the result depends on how to choose the Boltzmannized flux operator. 

In this paper, we have first considered a solvable model, i.e., the free particle case, to demonstrate the dependence 
of the rate constant on k intuitively. To study it, we have introduced a new scaling function and investigated its 
properties in detail. As a non-trivial case, we have discussed the (rr)-rate constant at very high pressure. Because 
under such conditions the reaction rate is determined by a propagator at short time, the ratio of the rate constant 
in the partial-split form of the Boltzmannized flux operator to that in the Yamamoto approach can be described 
in terms of the scaling function for the free particle case. We have found that the rate constant for recombination 
reaction in the Yamamoto approach is larger than that in the case of k ^ 1. In particular, the (rr)-rate constant in 
the Yamamoto approach is about twice as large as that in the Miller approach at a ~ 10. 

In conclusion, the partial-split form of the Boltzmannized flux operator is certainly an economical and powerful tool 
to perform numerical calculations for thermal rate constants of realistic reactions. However, for the recombination 
reaction it may underestimate the rate constant compared with the result calculated by the linear response theory. 
The difference could be seen if the experiments for recombination reactions could be performed under the conditions 
of very high pressure. 
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FIG. 1: Sketch of a potential surface in one-dimensional reaction versus the reaction coordinate s. 
denoted by r. 
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FIG. 2: One-dimensional schematic diagram of the potential for A + BC — » AB + C reaction. The compound region (ABC) 
is bounded by the dividing points on reactant (r) and product (p) sides. 
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FIG. 3; Same as Fig|5| but for the free particle case. 
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FIG. 6: Ratio of the (rr)-rate constants in the free particle case. The sohd (dot-dashed) [dotted] curve is for k = 0(0.5) [0.75]. 
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FIG. 8: Ratio of the (rr)-rate constants at very high pressure. The solid (dot-dashed) [dotted] curve is for k = 0(0.5) [0.75]. 



